Excitation gap of a graphene channel with superconducting boundaries 
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We calculate the density of states of electron-hole excitations in a superconductor-normal-metal- 
superconductor (SNS) junction in graphene, in the long-junction regime that the superconducting 
gap Ao is much larger than the Thouless energy Et — hv/d (with v the carrier velocity in graphene 
and d the separation of the NS boundaries) . If the normal region is undoped, the excitation spectrum 
consists of neutral modes that propagate along the boundaries — transporting energy but no charge. 
These "Andreev modes" are a coherent superposition of electron states from the conduction band and 
hole states from the valence band, coupled by specular Andreev reflection at the superconductor. 
The lowest Andreev mode has an excitation gap of Eo = ^(vr — \(j)\)ET, with cj) £ (— 7r,7r) the 
superconducting phase difference. At high doping (Fermi energy /i ^ Et) the excitation gap 
vanishes oc Eo{Et / fJ-)'^ , and the usual gapless density of states of Andreev levels is recovered. We 
use our results to calculate the ^-dependence of the thermal conductance of the graphene channel. 

PACS numbers: 74.45.-(-c, 73.20.At, 73.23.Ad, 74.78.Na 



I. INTRODUCTION 

The two-dimensional layer of carbon atoms known as 
graphene is a gapless semiconductor. A gap between con- 
duction and valence bands opens up if the layer is con- 
fined to a narrow channeli For a channel of width d the 
band gap 2Eo is set by the (ballistic) Thouless energy 
Et — hv/d, with v the (energy independent) velocity of 
electron and hole excitations in graphene. The size of the 
gap depends on the crystallography of the channel edges. 
In particular, for edges in the armchair configuration one 
has^ 



Eo = aEr, 



(1) 



with a = if the channel is a multiple of three unit cells 
across or a = 7r/3 otherwise. 

The interface with a superconductor provides an al- 
together different way to confine the carriers. At en- 
ergies below the superconducting gap Aq, the electron 
and hole excitations in a superconductor-normal-metal- 
superconductor (SNS) junction are confined to the nor- 
mal region. In usual metals this confinement leads to 
bound states known as Andreev levelsA^ They consist 
of counterpropagating electrons and holes converted into 
each other by Andreev retro-reflection at the NS bound- 
aries (see Fig.^). Andreev levels carry an electrical cur- 
rent (a supercurrent) across the NS interfaces, but they 
are "quasi-localized" along the interfaces. More precisely, 
the group velocity of the Andreev levels along the NS in- 
terface is much smaller than the Fermi velocity, and weak 
disorder fully localizes them& 

As pointed out in Ref.'?, Andreev reflection in undoped 
graphene is specular reflection instead of retro-reflection 
(see Fig.^). The consequences were investigated in that 
paper and in Ref. for a single NS interface. Here we 
consider the consequences for an SNS junction. 

We find that the transition from retro-refiection to 
specular reflection is accompanied by a transition from 
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FIG. 1: Three types of states in an SNS channel in graphene. 
The solid and dashed lines show the classical trajectories of an 
electron (filled circle) and a hole (open circle), converted into 
each other upon Andreev reflection at the superconductor. 
The transition from a localized Andreev level (a) to a propa- 
gating Andreev mode (b) occurs when the excitation energy e 
becomes larger than the Fermi energy ^ in the normal region. 
These two types of states are both charge-neutral. Purely 
electronic states (c) exist near grazing incidence. While the 
states of type (a) and (b) are sensitive to the phase difference 
cj> of the two superconductors, the states of type (c) are not. 



quasi-localized Andreev levels to propagating modes 
("Andreev modes"), provided that Et ^ Aq. This is 
the long-junction regime. (The states remain localized 
in the opposite short-junction regime Et 3> Aq, consid- 
ered in Ref. ffj) The transition (governed by the ratio 
fi/ET of the Fermi energy /j, in N and the Thouless en- 
ergy) has a drastic effect on the density of states. While 
the excitation spectrum is gapless for fi ^ Et, a gap 
opens up for fi < Et. The excitation gap 



Ea = ^{tt - \(t>\)Ej 



(2) 



has the same form as the band gap ^ for conflnement by 
armchair edges — with the phase difference S (— tt, tt) 
of the two superconductors taking over from the crystal- 
lographic phase a. 
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The Andreev modes have the same dispersion relation 
as the "armchair modes" for confinement by armchair 
edges, and they are also constructed out of states taken 
from two different valleys in the Brillouin zone. However, 
while the armchair modes contain cither electron states 
from the conduction band or hole states from the valence 
band, the Andreev modes are a superposition of conduc- 
tion and valence band states. As a consequence, the An- 
dreev modes transport energy but no charge along the NS 
interface — so they will play a role in thermal conduction 
along the interface but not in electrical conduction. 

The outline of this paper is as follows. The modes 
propagating along the channel are characterized by their 
dispersion relation in Sec. ^ Both exact numerical and 
approximate (but highly accurate) analytical results are 
given. From the dispersion relation we determine the 
excitation gap in Sec. lIIII and the density of states in Sec. 
IIVI contrasting in particular the low- and high-doping 
regimes. We derive the result (j2l for the excitation gap in 
the low-doping regime and show numerically that the gap 
closes oc E{){Et I IJ^Y' with increasing doping. One way to 
measure the gap is by tunneling spectroscopy. Another 
way, which we analyze in some detail in Sec. |V1 is by 
means of the thermal conductance of the channel (for 
heat flow parallel to the NS boundaries). We conclude in 
Sec. EH 



II. DISPERSION RELATION 
A. Quantization condition 

To calculate the dispersion relation of the An- 
dreev modes we solve the Dirac-Bogoliubov-De Gennes 
(DBdG) equation^ for the pair potential 



The three angles a± ,0±,j3 are functions of e and q, 



A(r) = 



Aq exp{i(j)/2) 


Aq exp(— i(?!</2) 



if a; < ~d/2, 

if -d/2<x< d/2, 

if a; > d/2. 



(3) 



We seek plane wave solutions \E'(a;,?/) = ip{x)e'^''''^ , with 
q the component of the wave vector parallel to the NS 
interfaces at x = ztd/2. The excitation energy e > of 
the mode is measured relative to the Fermi energy /i in 
the normal region \x\ < d/2. (The superconducting re- 
gions are assumed to be heavily doped, with Fermi energy 

The dispersion relation follows from the quantization 
condition derived from the DBdG equation in Ref. |^ 
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(5) 
(6) 



The quantization condition is invariant under /i - 
so without loss of generality we may take /i > 0. 

While Ref. dealt with the short-junction regime 
Et 3> Aq, here we are concerned with the long-junction 
regime Et <C Aq. (Since the Thouless energy Et = 
hv/d, the latter criterion is equivalent to the require- 
ment that the separation d of the NS interfaces is large 
compared to the superconducting coherence length ^ = 
Hv/Aq.) We furthermore restrict ourselves to low- lying 
excitations, e <C Aq. The relative magnitude of e and Et 
is arbitrary. For ease of notation we will use units such 
that hv = 1 in the intermediate calculations, restoring 
the units in the final results. 

For low- lying excitations e <C Ag the quantization con- 
dition simplifies to 



COS( 



cos 6* I cos 9_ + r sin 9+ sin 9_ = 0, 



where we have abbreviated 

1 -I- sin a+ sin a_ 



cos a+ cos Q!- 



(7) 



(8) 



The solutions to this equation can be represented in the 
form e = e^{q), where n = 0, 1,2, . . . is the mode in- 
dex due to the quantization of the motion in the x- 
direction and the superscript ± accounts for the different 
^-dependence of the modes. 



B. Exact solution 

The quantization condition |(7|) can be solved numer- 
ically. Results are shown in Figs. |21 El and ^ Only 
positive q is plotted, because s^{—q) = s^{q). 

The dispersion relation has three distinctly different 
branches, indicated in Fig. 2) corresponding to the three 
types of trajectories of Fig. 



(4) 



(a) The branch with hv\q\ < /i — e (red) describes in- 
traband electron-hole states, corresponding to the 
Andreev modes of Fig. The dispersion relation 
for these modes has small oscillations as a function 
of q on the scale 1 /d, around a smooth convex curve 
(see Fig. 13. 

(b) The branch with hv\q\ < e — fx (blue) describes in- 
terband electron-hole states, corresponding to the 
Andreev modes of Fig. ^p. The dispersion relation 
is concave without oscillations. 

(c) The branch with ?iw I g I > |e — /i| (green) corresponds 
to the purely electronic states of Fig. The hole 
component of the wave function can not propagate 
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10 20 30 40 50 60 70 

FIG. 2: Dispersion relation of the SNS junction, calculated 
numerically from Eq. Q for three values of the supercon- 
ducting phase difference <j) at h/Et = 0.1. The lowest modes 
en{q) with n = 0, 1, 2 are nearly degenerate for (j) = and 
nondegenerate for = -k/2 (thicker lines correspond to ej). 
For = TT all modes are nearly degenerate except the lowest 
one el ■ 



along the channel because the reflection angle q;_ 
of the hole is imaginary on this branch. The dis- 
persion relation on branch (c) is concave, without 
oscillations, and joined to branch (a) or (b) by a 
cusp singularity. 

After these exact results we continue with an ap- 
proximate, but highly accurate, analytical solution of 
the quantization condition. We consider separately the 
electron- hole modes with hv\q\ < \e — fi\ and the electron 
modes with hv\q\ > |e — 
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FIG. 3: Same as in Fig. El for fi/Er = 100. The smoothed 
dispersion relation l|14|l is indicated by dashed lines. 



C. Electron- hole modes 

For \q\ < |e — /i| (setting again ?it; = 1) we define the 
transverse momentum p by the relation 

e+-9^ = 2p/Et. (9) 

The solution to this equation is given by 

e=pJl~^^- (10) 
y - 

The condition \q\ < \e — ^\ is equivalent to |g| < 
= qc- The momentum q^ is the location of 
the cusp in the dispersion relation, beyond which the 
hole component of the mode vanishes. 
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If we express e in terms of p with the help of Eq. (|10|l , 
we can write 



(^2 _ _p2)2 _ (-^^^)2 



2/i 



? 9 ■ 



(lla) 



(lib) 



This allows to recast the quantization condition Q as 



COS( 



/ 2 _ 2 \ 2 

" \ ) [cos0 + cos(2p/i;T)], (12) 

which defines the quantization of the transverse momen- 
tum p = p^. 



For \q\ ^ Qc the solution to Eq. H12|l is given by 



(13) 



with n = 0, 1, 2, . . . and 4> € (— tt, tt). As g approaches the 
cusp at Qc the first term in Eq. (|12|l causes the dispersion 
relation e{q) to oscillate rapidly around a smooth curve 
s{q). This smoothed dispersion relation is obtained by 
substitution of Eq. p3|l into Eq. (|l()|l , resulting in 



ttEt [ n + - ± — 
' 2 2tt 



•1 - 



^2 _ (ttEj 



\q\ < Qc 



2 ^ 27r 



.(14) 



The smoothed dispersion for the lowest modes is indi- 
cated in Fig. O by a dashed line. 

To determine also the rapid oscillations, we proceed as 
follows. The quantization condition for p can be written 

as 

p^=nET(n+l±^), (15) 



27r, 

where the phases 7^ can be determined by iteration from 
Eq. (|12|l . The first iteration turns out to be already 
highly accurate in the high-doping regime 3> Et- It is 
given by 

In - 



(^2_p2)2COM 




(16) 



where the momentum p on the right-hand-side is taken in 
the zeroth approximation IjlSI) . The difference between 
the approximate analytical results of Eqs. HlUll5llt)|) and 
the exact numerical results plotted in Fig. |31 are not vis- 
ible on the scale of that figure. 




0.5 



FIG. 4: Dispersion relation of the SNS junction, calculated 
numerically from Eq. for <^ = and fj-fEr = 100. The 
curves show £t{q) with n = 5, 10, 15, . . . 60. The three types 
of states from Fig. are color coded; red=type a, blue=type 
b, green=type c. 



D. Electron modes 

For \q\ > |e — /i| the angle 9- becomes strictly imagi- 
nary. In this interval we define the transverse momentum 
p by 



2p/Ej 



(17) 



The condition \q\ > |e — /i| is then still equivalent to 
l^l > l/i — p2/yi| = From Eq. (|17|l we cast the branch 
\q\ > qc of the dispersion relation in the form 



q- 



4p2 



(18) 



where the momentum p = p^ is quantized. The exact 
quantization condition follows directly from Eq. (|7|l. 

For large longitudinal momenta \q\ 3> max {Et, qc) the 
reflection angle 9^ w zLiq/Ex of the hole takes on large 
imaginary values. Therefore both sin 9- and cos 9_ in the 
quantization condition (0 are exponentially large and 
the (/)-dependence of the solution can be neglected. This 
shows that the electron modes are insensitive to the su- 
perconducting phase difference across the channel. 

At the cusp \q\ — qc of the disperrsion relation we find 
6'_ = and e = The coefficient r in Eq. tends 

to infinity, leading to 



lim r sin ( 

kl->9c 



pEt 



(19) 



The quantization condition at the cusp thus simplifies to 



cos <p + cos(2p/ Et) ^ ^ ( - - -] sm(2p/ET). (20) 

For /i ^ Et the condition is reduced to sm{2p/ET) — 
with solution 

p+ = nET{n + 1), p- = ttEt {n+^), (21) 
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FIG. 5: Double-logarithmic plot of the energy dependence of 
the excitation gap, calculated numerically from Eq. (|7|l, for 
three different values of the superconducting phase difference 
cj> (colored lines). The straight black line is the asymptote 



FIG. 6: Thermodynamic limit Et — hv/d ^ of the density 
of states p of the SNS junction, according to Eq. (I24I I with 
po ~ 4:£/nhv. The total density of states p (solid) is the sum 
of the density oc /i + e of electron states (dashed) and the 
density cx |p — e| of hole states (dotted). 



which is, again, independent. The quantization con- 
dition formally corresponds to 7+ = tt, 7" = in 
Eq. (O. 

III. EXCITATION GAP 

At small doping pL ^ Et the excitation gap is given by 

which is the energy of the lowest mode at q = 0. For 
/i ^ Et the gap is given by the energy of the lowest mode 
at a nonzero longitudinal momentum \q\ ^ Qc, which cor- 
responds to the deepest minimum of the oscillatory dis- 
persion relation. We have not succeeded in determining 
this minimum analytically from the quantization condi- 
tion Q, but we have a very accurate numerical solution. 

Results for different values of </> are presented in Fig.O 
One can see that the ratio £'gap/i?o('/') depends only 
weakly on the superconducting phase difference (f> and 
that the crossover to a decay oc happens in a narrow 
interval around pi = ttEt- As shown by the black line in 
Fig. |S1 the large-/x asymptote is given by 

i?gap=c(0)i?o(0)(7ri?T/M)', (23) 
with c((/)) increasing from 1/2 at = to 1 at |0| = tt. 

IV. DENSITY OF STATES 

A. Thermodynamic limit 

Before turning to the calculation of the density of 
states at finite Et , it is instructive to first determine the 



behavior in the thermodynamic limit d 00 ■i^ Et — > 0. 
The DBdG equation then decouples into separate Dirac 
equations for electrons and holes. The total density of 
states p{e) — Pe(e) + Phi^) is the sum of the electron 
density of states pds) oc p + s and the hole density of 
states Ph{£) oc |/i — e|, leading to 

ACd 

P(^) = max(^, e), if e > Et- (24) 

Here C is the extension of the junction in the y-direction 
and the factor of 4 accounts for the spin and valley de- 
generacies. (For a derivation of Eq. 1)24(1 directly from 
the quantization condition lO, see Add. 1X1) 

In Fig. |Hlwe plot the density of states p of the DBdG 
equation together with the separate electron and hole 
contributions pe and ph- The superconducting prox- 
imity effect will introduce fine structure in p on the 
scale of the Thouless energy Et, as we will determine 
in the next subsections. We consider separately the low- 
doping regime p ^ Et, where the contribution from in- 
ferband electron- hole modes dominates, and the high- 
doping regime p ^ Et, where the intrahand electron- 
hole modes dominate. 



B. Low-doping regime 

To determine the excitation spectrum in the low- 
doping regime, we take the p ^ limit of Eq. I|14|l . 
resulting in 

e± = ^J{nvqy + iTTET)Hn+^±cl>/27Ty, (25) 

for n = 0, 1, 2, . . . and G (— tt, tt). (There is no need to 
distinguish e from e, because the dispersion relation does 
not oscillate in this regime.) The two series of modes 
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FIG. 7: Density of states of the SNS junction in the low- 
doping regime, for superconducting phase difference = 
(sohd curves) and (f> = -n (dashed curves). The curves are 
calculated from Eq. (1261 . normalized by po = 'iC/nTiv. The 
excitation gap for = is at i5o = -eEtII- 



FIG. 8: Smoothed density of states of the SNS junction in 
the high-doping regime for </> = 0, calculated from Eq. 13111 . 



C. High-doping regime 



and e„ are nondegenerate, except for i 
lowest mode £q is nondegenerate also for 
In Fig. [3 we plot the density of states 



0,7r. (The 



Pie) 



n=0 ± 



dq 



■kTivEt 



5:5:(x±)-v2e(x±), (26) 



=0 ± 



(27) 



with Q the unit step function. The excitation spectrum 
has a gap at the energy Ea given by Eq. (jSJ. The gap 
closes for |(^| = tt, when p = AC/nhv = po is constant 
at low energies. At large excitation energies e ^ Et the 
sum over n in Eq. 1)2611 may be replaced by an integral, 
resulting in a linearly increasing density of states. 



Pie) 



4Ce 
ttHvEt ' 



if e > Et, 



(28) 



in agreement with the thermodynamic limit (|24(l . 

The group velocity in the y-direction of the n-th 
mode is given by the derivative 



hdq 



(29) 



For each propagating mode f ^ — > w with increasing ex- 
citation energy. These are all interband electron-hole 
modes. The purely electronic modes are pushed to 
\q\ E^/p — > oo in the low-doping regime p/E^ — > 0, 
while the intraband electron-hole modes can not propa- 
gate if e > p. 



As shown in Sec. Ill CI the electron-hole branch of the 
dispersion relation for p 3> Et is a rapidly oscillating 
function. Small local variations in the separation d of 
the NS interfaces, on the scale of the Fermi wave length 
Xp = hv/p, will average out these oscillations, leaving 
the smoothed dispersion relation H14|l . In the large- /j, 
limit this reduces to 



ttEt 
P 



^±±.]^p2-(nvq)^, (30) 



The branch of purely electronic states (for hv\q\ > p) is 
not described by Eq. (jSOJ, but since it contributes neg- 
ligibly to the density of states for e p we need not 
consider it here. 

The smoothed density of states is given by 



Pie) 



AC 

TT 



EE 

n=a ± 
ACpi 



dq 



- ,.,„sjEi:(^-.f)-"^e(r.f) 

-< 71=0 ± 

X (n+i±0/27r)-\ 
r± = 7T\n+^±<j>/27r)'-{e/ET)'. 



(31) 
(32) 



We plot it in Fig. for (/) = 0. 

The peaks in the density of states at — irETin + 
i ± 4'/2tt) are analogous to the De Gennes-Saint James 
resonancesiS in conventional Josephson junctions. The 
lowest resonance is at the same energy Eq = ^(Tr— \(t>\)ET 
as the gap Q in the low-doping regime — however, in 
the high-doping regime the density of states is gapless, 
vanishing linearly at small excitation energies with a 0- 
dependent slope: 



Pie) 
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nhvEl cos2(0/2)^ 



if £ < Eq. 



(33) 
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FIG. 9: A temperature diflterence AT = Tl - Tu between 
the two ends of the graphene channel drives a heat current 
Iq, carried by Andreev modes in the normal region at tem- 
peratures below the gap Aq in the superconductors. 



The slope diverges when — > tt, because then the lowest 
resonance is at e = 0. At high excitation energies e ^ Et 
(but still e fi) the density of states approaches a 0- 
independent limit, 



■kTivEt ' 



if £ > Et, 



(34) 



(EH). 

* — de^/hdq corresponding to 



in agreement with Eq 

The group velocity v, 
the smoothed density of states is of order vEt/^ 
vXp/d for hv\q\ <^ fJ., much smaller than the carrier ve- 
locity V. This as expected from the classical trajectories 
in Fig.HJi. 



V. THERMAL CONDUCTANCE 

The thermal conductance Gthermai = Iq/^T of the 
graphene channel, for heat flow Iq parallel to the NS 
boundaries, can be measured by applying a tempera- 
ture difference AT = Tl — Tr between the two ends 
of the channel (see Fig. O. Experiments of this type 
have been performed in metals by Eom, Chien, and 
Chandrasekhar^^ and analyzed theoretically in Refs. 0, 

ES 

To determine the thermal conductance of the graphene 
channel we start from the Landauer-type formula— 
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FIG. 10: Thermal conductance of the SNS junction in the 
low-doping, low-temperature regime {fi,kBTo <^ Et), cal- 
culated as a function of the superconducting phase differ- 
ence (j> (modulo 2tt) from Eq. 13711 . The peak value equals 
Gpcak — Tvk%To/3h. 



In the low-doping, low-temperature limit ^, ksTo ^ 
Et only the lowest mode contributes and the thermal 
conductance reduces to 



G 



thermal 



k%To 



dx ■ 



Eo/ksTo cosh (x/2)' 



(37) 



with Eq the gap given by Eq. As illustrated in Fig. 
1101 the thermal conductance in this low-doping, low- 
temperature regime vanishes unless the superconduct- 
ing phase difference is in a narrow interval of order 
ksTo/ET around tt (modulo 2tt). The peak at (j) — n has 
height 



G 



peak 



3 h 



(38) 



G 



thermal " 



/ rfe£^§^Xl^" ('^)' if we replace £;o by 

•'O n low temnerature resi 



In the high-doping limit ^ 3> Et we may distinguish 
a moderately-low temperature regime E^/^ <C fcsTo ^ 
Et and an ultralow temperature regime fcsTo ^ E^/^. 
In the ultralow temperature regime it is again only the 
lowest mode which contributes, so Eq. H37|) remains valid 



from Sec. mil In the moderately- 
low temperature regime there remains a large number 



valid for small temperature differences AT <ti TttTr- 
(The factor of 4 is again from the spin and valley de- 
generacy.) We assume that the mean temperature Tq = 
{Tl + Tif)/2 is much less than Ao/ks, so that the ther- 
mal current through the superconductors is suppressed 
exponentially.^ The function /(e) = [1 -I- exp(e/fcBTo)]^^ 
is the Fermi function and T„ is the transmission proba- 
bility of the n-th propagating mode along the channel. In 
a ballistic channel each of the N{e) propagating modes 
at energy e has transmission probability 7^ = 1, so we 
obtain the thermal conductance 



N{e) 



ttEt 



(39) 



of modes that contributes at energies e ^ ksTo. Substi- 
tution into Eq. H36(l gives the thermal conductance 



G. 



thermal 



2.34 



I fj.kBTo 
El " 



(40) 



G 



1 



thermal 



27r?ifcsT(f 



de ■ 



coali {e/ 2k bTq) 



(36) 



The thermal conductance is insensitive to the supercon- 
ducting phase difference because of the vanishing excita- 
tion gap in the high-doping regime. 
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VI. CONCLUSION 



compactly written as 



We have shown that a graphcne channel with supercon- 
ducting boundaries supports a type of propagating modes 
along the channel that do not exist in usual SNS junc- 
tions. These "Andreev modes" exist because the Andreev 
reflection close to the Dirac point of vanishing Fermi en- 
ergy fi is specular^ The Andreev modes are charge neu- 
tral, so they transport energy but no charge along the 
channel. 

The thermal conductance due to the Andreev modes 
depends strongly on the superconducting phase differ- 
ence (j), because of the 0-dependent excitation gap Eq 
of the Andreev modes. Away from the Dirac point the 
character of the Andreev reflection changes from specu- 
lar reflection to retroreflection. The excitation gap closes 
and the thermal conductance becomes 0-independent. 

The closing of the excitation gap with increasing dop- 
ing can be studied directly by point contact spectroscopy 
(tunneling into the graphene layer via a tunnel probe on 
top of the layer). 
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4(9) 



— 



q < Im-pVmI, 



(Al) 



V^^Tlf-H, q>\fi-p^/fi\, 



where the momentum p = is quantized according to 
Eqs. (|15|l and (|21|) . The density of states is given by 



n=0 ± 



dqd{e - etiq)). 



(A2) 



Since the quantization condition of the momentum p is 
linear in the mode index n, we can replace the summation 
over n with the integration over p in the thermodynamic 
limit d — > oo. In this limit we can ignore the dependence 
of the phases 7^ on q. The integral Eq. (jA2|l results in 



Pie) 



8Cd 



dp 



8£,d 
4Cd 



min|e,^/e7l} 



e hi^ — p^ 
P\ p^ — 



dp 



v/(e + /i)2 - V 



max (e, /x). 



(A3) 



(A4) 



APPENDIX A: DERIVATION OF EQ. FROM 
THE QUANTIZATION CONDITION 

It is instructive to calculate the density of states in 
the thermodynamic limit directly from the quantization 
condition ((TJ. The dispersion relation for all q can be 



The first and second integral on the right-hand-side 
of Eq. (|A3p are, respectively, the contributions from 
electron-hole and electron modes to the density of states 
in the thermodynamic limit. Even though each integral 
is a non-trivial function of energy, their sum reduces to 
the elementary result (|24|l . confirming the consistency of 
our analysis. 
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